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Abstract 


The  synergistic  effects  of  an  adspecies/surface 
system  influenced  by  a  coherent  laser  field  and  an  in¬ 
coherent  thermal  phonon  field  are  theoretically  inves¬ 
tigated.  A  generalized  master  equation  is  derived  in 
the  Schrodinger-Markoff  picture  and  the  total  stimu¬ 
lated  transition  rates  are  decomposed  into  three  parts: 
thermal  phonon,  laser  and  the  interference  terms.  The 
energy  transfer  dynamics  is  pictured  via  the  evolution 
equation  of  the  average  excitation  of  the  active  mode. 

A  random  phase  of  the  off-diagonal  matrix  elements  in¬ 
duced  by  the  laser-stimulated  surface  processes  is 
introduced.  Finally,  possible  applications  of  the 
laser/phonon-stimulated  master  equation  are  discussed. 
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I .  Introduction 


Laser-matter  interactions  of  heterogeneous  systems  have  been 
recently  studied  and  the  possible  applications/implications  of  laser- 
stimulated  transient  effects,  such  as  materials  processing  and 
evaluation,  semiconductor  lithography  and  microelectronics  and  the 
development  of  new  lasers,  motivate  the  fundamental  studies  in  this 
newly-developed  field. 

Laser-stimulated  surfaces  processes  (LSSP)  such  as  migration,  de¬ 
sorption,  adsorption,  ionization,  recombination  and  predissociation 
usually  combine  the  laser  transient  effects  and  the  subsequent  energy 
relaxation  of  the  excited  adspecies . ^  Due  to  the  multiple  time  scales 
of  the  heterogeneous  system,  LSSP  may  be  grouped  into  several  types, 
e.g.,  mode-selective,  adspecies-selective,  migrational-selective  and 
local-heating,  etc.  These  selective  and  nonselective  aspects  are 
characterized  not  only  by  the  coherent  properties  of  the  laser  rad¬ 
iation,  frequency,  polarization,  intensity  and  duration,  but  also  by 
the  pumping  rate  compared  to  the  relaxation  or  coupling  rates  among 

inter-  and  intra-modes,  e.g.,  electron-phonon  interaction,  phonon- 

2 

phonon  coupling  and  active-mode-phonon  coupling. 

In  a  previous  paper, ^  we  studied  the  desorption  dynamics  by  means 
of  a  laser-stimulated  rate  equation,  where  t>e  thermal  effects  and 
the  laser  excitation  were  assumed  to  be  two  independent  driving  forces, 
and  we  focused  on  a  cold  surface  where  the  laser  excitation  plays  the 
dominant  role.  In  the  present  paper,  we  shall  study  the  synergistic 
effects  of  laser/phonon  excitation  where  both  the  laser  and  thermal 
phonons  may  play  equally  important  roles. 
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II .  Equation  of  Motion  in  the  Schrodinger  -  Markoff  Picture 

We  consider  a  model  system  consisting  of  a  group  of  species 
adsorbed  on  a  solid  surface  and  subject  to  laser  radiation.  The 
total  vibrational  Hamiltonian  of  the  system  may  be  expressed  in  the 
following  form: 

8  ■  HA  +  Hb  +  +  Hbb  +  Hab  +  +  Hbf.  (1) 

and  Hb  are  the  unperturbed  Hamiltonian  of  the  optically-active 
(A)  mode  and  the  phonon  bath  (B)  modes  respectively;  H  and  H__ 

AA  0X3 

represent  the  intramode  couplings  of  the  active  and  phonon  bath  modes, 
respectively;  H^  and  HBp  represent  the  A  and  B  mode  coupling, 

the  A  mode  and  laser  field  and  the  B  mode  and  laser  field  interactions, 
respectively.  For  the  situation  where  the  vibrational  spectrum  is 
partitioned  into  a  high-frequency  region  consisting  of  the  A  modes 
resonantly  excited  by  a  laser,  we  may  neglect  Hgp  since  the  low-frequency 
B  modes  are  far  off-resonant  with  respect  to  the  laser  radiation.  A 
second-quantization  representation  of  the  total  Hamiltonian  was  pre¬ 
viously  shown  in  Ref.  3  consisting  of  the  combinations  of  the  ladder 
+  + 

operators  a  (a)  for  the  active  mode  and  b  (b)  for  the  phonon  modes. 

For  example,  H  =  V{t)  (a  +  a)  where  V(t)  is  proportional  to  the 

AT 

electric  field  of  the  radiation. 

In  the  Schrodinger  -  Markoff  picture  (SMP) ,  the  equation  of 
motion  for  the  reduced  density  operator  S  =  Trfip,  traced  over  the 

phonon  modes,  with  the  total  Hamiltonian  given  by  eg.  (1)  is  found 
3,4 


to  be 


||  =-i[0)effata,s]  -  iV(t)([a+,S]  +[a,S]) 


ef  f 


+  ^(2aSa+-a+aS-Sa+aj  +  y^j  (a+Sa+aSa+-a+aS-Saa+) ,  (2. a) 


with 


%££<« 


=  o)Q  -  e*(aTa+l)  +  ^(t). 


(2.b) 


y,  *  2 it  E|K 


v,  nv  6(w0-2e*n-  Z  u>j)  , 


(2.c) 


n  =  n(H.+l)  -  II  n .  , 
v  t  j  -1  J 


(2.d) 


3  J  3 

where  (ujg-e*)  and  are  the  fundamental  frequency  of  the  active  mode 

and  the  j-th  phonon  mode  with  the  vibrational  quantum  number  n  and  Rj, 

respectively.  The  important  features  of  the  above  equation  of  motion 

in  SMP  are:  (i)  the  many-body  effects  of  the  thermal  phonon  modes  are 

effectively  reduced  to  a  multiphonon  dampinq  factor  y^  which  is  char- 
* 

acterized  by  the  mean  occupation  number  n^  (at  temperature  T) ,  the 
coupling  strength  Kv  and  the  delta  function  5(u0,u)j)  for  conserving 
energy;  (ii)  the  active  mode  is  treated  as  a  stochastic  anharmonic 
quantum  oscillator  with  an  effective  frequency  u>eff  given  by  an 
anharmonic  correction  e*(a+a+l)  and  a  random  phase  or  frequency  shift 
<j>(t);5  (iii)  the  dephasing  effects  of  the  active  mode  may  be  caused  by 
a  laser/surface-induced  mechanism  such  as  thermal  fluctuation  of  the 
active  dipole  moment,  migration  or  dislocalization  of  the  adspecies, 
phase  difference  between  the  laser  (coherent)  and  phonon  (incoherent) 
fields  coupled  to  the  active  mode,  random  imperfection  of  the  solid 
surface  and  finally,  the  effects  of  the  intramode  couplings  and 


Ii 


BB* 
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The  evolution  equation  of  the  density  matrix  elements  for  the 

4  6 

vibrational  states  (m,n) ,  ' 


-ir  *  TK  l  (Hmk  skn  '  smk  Hkn>  '  W  (S  '  5>»n  ' 
can  be  written  for  our  system,  described  by  the  operator  equation 
[eq.(2)],  as 

3S 

=  -i  t(o)0-e*+<)>)(m-n)-e*  (m2-n2)-iYQ] 


(3) 


[  jm  S^n  -  fn+I  W.tJ5i  sm+l,n 

+  T  [  2  'J<m+1)  (n+1)SmH,n+l  *  (m+n)sm,n)] 

+  Vj  [/■"“  sm-t,n-l  +  ^m+1)<n+1>Sm+1,n+l  ' 


-F 

(",+n+1)Sm,n]' 


(4) 


where  y ’  is  a  quantity  proportional  to  the  derivative  of  the  dipole 
transition  matrix  element  between  ground  and  first-excited  states, 

E ( t)  =  Egcoswt  is  the  electric  field  of  the  linearly  polarized  laser, 
and  Y0=i(Ym+Yn)  is  the  mean  natural  width  for  the  radiative  decay  from 
state  m  to  state  n.  . 

Ill .  Laser/Phonon-Stimulated  Master  Equation 

The  evolution  equation  of  the  energy  population  is  given  by  the 

diagonal  term,  P  m  in  eq.  (4)  , 

n  nn 


dP 


-dr  =  [  (n+1)  (vl)  Pn+1  +  niij  pn-l  -(2nSj+n+Bj) 

-(iy'E(t^[^(Sn.1;n  -  Snin.1  -  JZZl  (Sn>n+1  - 


pn]-YnPn 

Vl,n>]  '(5) 
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where  pn±i"sn±^  n±1  and  the  off-diagonal  terms  Sn,  n±1  can  be 
calculated  iteratively  by  means  of  the  equation6 

3s<j> 

_BS m  -i.  ru  SI  -  v  t S— S )  (6) 

-Tr-  -  [H,  SJmn  ymn  (S  S)^  ,  (6) 

with'the  condition  (t)  ■  sj;?*  (t)S__  and  SjfO)  «0  for  m*n,  we  obtain 

mn  ran  ran  ran 

Sn^n+l(t)  -  -(fljn+l  J  u’E(f)exp  [j  ft  (t")  (t'-t")  df]  * 


tPn+i(t,)'Pn0)  (t’)]  dt'  » 


(7. a) 


where 


ftttJ-itwQ-Ze^vn+l)  +<}>(t)]  -  ( 2n+l-4iin+4n)  Yj/2  -y0 


(7.b) 


Note  that  ft(t)  contains  a  random  phase  factor  <fr(t) •  By  assuming  a 

"white  noise",  <$(t)>*0,  and  a  Markoff ian  correlation,  <<fr(t)<Mt')>  * 

2Y2S(t-t'),  we  may  evaluate  the  ensemble-averaged  quantity  [over  the 

random  variable  <J>(t)],  as 

t'  »  (-Y2t)X1  “Y2t 

<  exp  [i(  4>(t")  (f-t")dt"]>  *  Z0~HT -  *  e  ’  (8^ 


which  is  governed  by  the  dephasing  (T2)  factor  y2*  Equation  (7. a) 


then  takes  the  form 


:Sitn+l(t)>“^in+1  $  u'EU’Jexplft^  [<?£{  >  >“<Pn°)  (t')>]dt»  , 


(9. a) 


ft^(t-t')  ■  ^ i  [Wq-2e* (n+1 (t-t') 
F1  “  Y1(2n+l+4nn+4n) /2  +  Y2  +  Yq  • 


(9.b) 


(9  .c) 
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By  the  same  procedure#  we  find 

a  a  ® 

<Sn-l#n(t)>  SoP'E(t,)exp[S2l  (<P^°[(t')  >  -  <P^0)  (t  • )  >]  dt'  ,  (10. a) 

n2(t-t')  -  (i[to0-2e*n]-r2)  (t-t')  #  ,  (10. b) 

r2  “  Y1(2n-l+4nn)/2  +  Y2  +  Y0  *  KlO.c) 

Substituting  eqs.  (9)  and  (10)  into  eq.  (5),  we  obtain  the 
ensemble-averaged  generalized  maser  equation  (GME) #  neglecting  the 
ensemble-averaged  notation  <...>,  for  the  first-order  solution# 


dPn 

ST  "  Yi  [<n+1)  (5j+X)Pn+l  +  n5jpn+l  '  (2»5j+n+5j)Pn]  -YnPn 


iV'  (Wn,n-X(t-t,>tpn-l(t'>  -pnit,>l  +  [pn+l<t'>-pn(t’>  l)  , 


with  the  time  dependent  transition  rates 

W  (W) -  2  frTOlu'AfEMEM  (*-?>] , 

)-Z7iU<A.fRQE(fi  . 


The  above  GME  may  be  reduced  to  the  Pauli  master  equation  ,ME)  by  using 
the  fact  that  the  population  inversion  (Pn±1(t) -Pn(t) ]  is  usually  a 
slowly  varying  function  in  the  time  of  and  hence  can  be  factored 

outside  the  integral  in  eq.  (11) .  Taking  an  average  over  the  cycle  of 
of  the  field,  2tt/w  ,  and  assuming  a  near  resonance  condition  such  that 
the  integrated  transition  rates  W  . (t)  become  time-independent,  and 
we  obtain  a  simple  laser/phonon- stimulated  master  equation  (LPME) 
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dt  (Wn,n+l+Wn^n+l)(Pn+l"Pn)  +  (wn,n-l+wn,n-l)  (Pn-l“Pn) 

+  (Rn,n+lPn+l  “  Rn,n-lPn)  “YnPn  , 
where  the  stimulated  transition  rates  are  given  by 

^.n+l  ■  rl5j(n+1)-  , 


w* 

n,n-l 

n,n+l 


n,n-l 

g1>2(S) 


finjn 


5(^)2 .,«!»  . 

-  r1#2/<^l,2+ri,2)  ' 


“  w0“2e* (n+1) -u>  , 
22  -  o)Q-2e*n-w  . 


(13) 


(14. a) 
(14. b) 
(14. c) 

(14. d) 

(14. e) 
(14. f) 
(14. g) 


The  phonon- induced  relaxation  rates  R  ...■(n+l)Y1  and  R  „  .•my, 

n,n+i  l  n,n~i  i 

describe  the  irreversible  energy  dissipative  processes  for  excited- 
species  on  a  cold  surface.  The  final  terra  in  eq. (13)  describes  the 
spontaneous  radiative  decay  of  the  n-th  vibrational  levels. 


IV.  Discussion  and  Conclusion 

It  is  shown  in  the  LPME  [eq.(lO)]  that  the  active  mode  is  in¬ 
fluenced  by  a  coherent  laser  field  and  an  incoherent  thermal  phonon 

field  with  total  stimulated  transition  rates  (vP  +  wf^^..)  which  may 

n,m±  n,  nil 

be  further  decomposed  into  three  components  (V^  ^ ,  +W1,  ,  -W1  , )  , 

n,ntl  n,n±l  n,ntl 

corresponding  to  phonon-stimulated,  laser-stimulated  and  the  "inter¬ 
ference"  transition  rates,  respectively.  wj^n+1  is  given  in  eq.  (14) 
and 


n,n+l 


(15. a) 


1 
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n,n-l 


(15. b) 


y< 


+  Yo>, 


*1,2  '0/  '“1,2 

are  simply  the  surface-free  laser  pumping  rates  characterized  by  the 
Rabi  frequency  and  a  Lorentzian  with  FWHM  equal  to  the  natural  level 
width  Yq*  The  interference  term  is  related  to  W*1  and  W1*  by 


(15. c) 


<n±l  =  ’ 


(16) 


for  Yq  >>  Ni  2  Y1  +  Y2'  ^1  2  * 

For  a  monatomic  lattice  with  phonon  frequency  «  Debye  frequency 

(wD) ,  the  multiphonon  coupling  factor  [eq.  (2.c)j  may  be  approximated 
7  N—  1 

as  Yj_  s  Y^(0)T  [exp(ho)0  /kT)  -1]  ,  with  the  zero  temperature  factor 
Y1(0)  ~  a  ,  where  N«  (a)Q-2e*n) /wD  in  the  order  parameter  and  0<a<<l. 
For  a  highly  excited  active  mode,  the  multiphonon  order  parameter,  N, 
decreases  rj  the  anharmonic  correction  2e*n  increases  and  in  turn 
results  a  stronger  coupling  factor  y^*  A  stronger  coupling  factor 
Y^  results  in  a  significant  interference  term  W1  as  shown  in  eqs. 

(16) .  Therefore,  the  total  transition  rate  in  the  LPME  can  be 
approximated  by  the  sum  of  W  and  W  only  when  W  is  negligible, 

i.e.,  2y1+y2^  <<y0  . 

To  give  a  clearer  picture  of  the  synergistic  effects  of  laser/ 

phonon  excitation  processes,  we  consider  the  average  excitation  of  the 

t  t 

active  mode,  <nA>*<<a  a>>  =  <Tr(Sa  a)>,  where  <<•••  >>  denotes  an 
ensemble  average  over  the  A  and  B  modes.  The  evolution  of  <nA>  can 
be  obtained  from  eq.  (2)  as 

<nA>  =fv«#E(tyk]  [<<a  >>  -  <<a*>>]  +  -  Y1<nft>  ,  (17) 


which  shows  that  the  active  mode  is  pumped  by  two  forces,  a  laser 
field  u'E(t)  and  a  thermal  phonon  field  y^n^,  while  it  is  damped  by 
a  multiphonon  relaxation  rate  the  last  term  in  eq.  (17) .  For 
the  cage  of  no  laser  field,  E(t)-0,  <nA>  reaches  the  steady-state 
value  rL  and  is  in  thermal  equilibrium  with  all  the  other  modes. 

With  laser  radiation,  a  nonequilibrium  state  where  the  active  mode 
has  a  higher  vibrational  energy  than  the  others  is  provided,  and  se¬ 
lective  excitation  is  possible  in  a  time  scale  shorter  than  that 
of  the  multiphonon  relaxation 

In  conclusion,  we  propose  some  possible  applications  of  LPME 
[eq.  (13) ]  where  the  laser  and  the  surface  play  equally  important 
roles:  (i)  solving  the  LPME  for  Pr  by,  e.g.,  generating  function 
method,  enables  us  to  calculate  physical  quantities  such  as  the  de¬ 
sorption  probability1"  and  the  desorption  rate  (or  the  mean  first 

Q 

passage  time) ;  (ii)  the  adspecies/surface  spectrum  governed  by 

the*  correlation  function  of  the  off-diagonal  matrix  elements  [eqs. 

(9)  and  (10)]  will  provide  us  surface  dynamical  information  such  as 

the  surface  coverage  and  the  migration  rate  of  the  adspecies  through 

the  dephasing  broadening  of  the  spectrum  Y2  which  reflects  the  col- 

lisional  effects;  (iii)  replacing  the  vibrational  population  function 

by  the  lattice  site  probability,  the  LPME  in  eq.  (13)  then  may  be 

used  to  describe  the  dynamical  motion  of  the  adspecies  on  the  solid 

surface  in  which  the  vibrational  pumping  rates  may  be  regarded  as 

9 

the  lattice  hopping  rates. 
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